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Clustering and gelation of hard spheres induced by the Pickering effect 

Andrea Fortino 

Theoretische Physik II, Physikalisches Institut, Universitat Bayreuth, 

Universitdtsstrafle 30, D-9544'^ Bayreuth, Germany 

A mixture of hard-sphere particles and model emulsion droplets is studied with a Brownian 
dynamics simulation. We find that the addition of non-wetting emulsion droplets to a suspension 
of pure hard spheres can lead to both gas-liquid and fiuid-solid phase separations. Furthermore, we 
find a stable fluid of hard-sphere clusters. The stability is due to the saturation of the attraction 
that occurs when the surface of the droplets is completely covered with colloidal particles. At 
larger emulsion droplet densities a percolation transition is observed. The resulting networks of 
colloidal particles show dynamical and mechanical properties typical of a colloidal gel. The results 
of the model are in good qualitative agreement with recent experimental findings [Koos, E. and 
Willenbacher, N. (20f f ). Science, 33f , p. 897] in a mixture of colloidal particles and two immiscible 
fluids. 

PACS numbers: 61.20.Ja,83.80.Hj,81.16.Dn 



Controlling the flow and the mechanical properties of 
complex fluids has important implications in both science 
and technology. Therefore, the possibility to change the 
viscoelastic and flow properties of colloidal or nanoparti- 
cle suspensions simply by adding a small amount water 
is intriguing in its simplicity and elegance. This sur- 
prising result was obtained experimentally by Koos and 
Willenbacher [l| (see also Ref. H). They found that the 
addition of an immiscible non-wetting fluid to a stable 
colloidal suspension leads to the formation of emulsion 
droplets that act as a bonding agent between the col- 
loidal particles. As observed in microscopy images, the 
transition from a fluid to a gel is driven by the formation 
of a network of colloidal particles bonded by emulsion 
droplet. It is well known that colloids trapped at droplet 
interfaces reduce the interfacial free energy (the Picker- 
ing effect [3] ) , and this basic physical mechanism is used 
to produce colloidal clusters Q-Q as well as to stabilize 
emulsions [7|, Q and bigels Q • 

A detailed understanding of the gel-transition mecha- 
nism driven by a low density of secondary fluid emulsion 
droplets is missing at present. The aim of this article is 
to introduce a minimal model that describes the network 
formation when non-wetting emulsion droplets are mixed 
with a stable suspension of purely repulsive colloidal par- 
ticles. 

We find that the Pickering effect drives the cluster- 
cluster aggregation that is behind a transition from 
a fluid of clusters to a percolated network, and that 
the properties of these networks are those of colloidal 
gels |10l - KL3| . Since the Pickering mechanism saturates 
once the droplet surfaces are completely covered with col- 
loids, a fluid of stable clusters is achievable without the 
addition of a long range repulsion [lj,[ia]. In addition we 
find that the Pickering effect can induce gas-liquid and 
gas-solid phase separations. 



We describe the suspension as a binary mixture of Nc 
colloidal particles with hard-sphere diameter CTc and Ad 
spherical droplets of diameter crd- The total interaction 
energy U is the sum of colloid-colloid, droplet-droplet 
and colloid-droplet interactions. 
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Nri 



= ^0cc(|rj -rj|)-t-^<?!)dd(|Rj-Ri|) 
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+ ^^0cd(|r. -Rj 



(1) 



where k^ is the Boltzmann constant, T is the tempera- 
ture, r^ is the centcr-of-mass position of colloid i, Rj is 
the center-of-mass position of droplet j, (j)cc is the colloid- 
colloid pair interaction, (jjcd is the colloid-droplet pair in- 
teraction, and (j)dd is the droplet-droplet pair interaction. 
The colloid-colloid and droplet-droplet interactions are 
hard-sphere like 
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otherwise, 
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where i ^ d, for the droplet-droplet interaction and i = c 
for colloid-colloid interaction. The steep repulsive con- 
tinuos potential was chosen to allow integration with a 
Brownian Dynamics algorithm |16l |. The repulsion be- 
tween droplets is aimed at modeling the repulsion of 
charged droplets; with this assumption droplets coales- 
cence can also be neglected. Furthermore, the shape of 
the droplets remains spherical. 

The colloid-droplet interaction is aimed at modeling 
the Pickering effect. The loss of interfacial energy Q 
when a particle is trapped at the surface of a droplet is 
modeled simply by a parabolic well of depth e 
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where a=(— e+ ^Jt 
B = i<7d + a,)/2- 
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Figure 1. (Color online) Equilibrium phases at (fjc = 0.1. 
Squares indicate homogeneous fluid phases, circles indicate 
gas-liquid coexistence, while triangles indicate fluid-binary 
crystal coexistence, a) Droplet size o-d/o"c=0.75. b)Droplet 
size (7d/crc=1.5. 



contact angle between droplets and colloids to 6 = 150°, 
and 138°, for crd/cc^l-S, and 0.75, respectively. Cur- 
vature and merging of droplets are neglected. A more 
detailed model Q would be necessary to describe all pos- 
sible droplets behaviors. 

We carried out Brownian Dynamics (BD) simula- 
tions [16| (hydrodynamics interactions are neglected) 
with time stcptJt = 1 x 10~^tb, with re = a^/Do, where 
Dq is the Stokcs-Einstein diffusion coefficient. The ran- 
dom forces mimic the interaction between particles and 
solvent, and are sampled from a Gaussian distribution 
with variance iD^St. The system is at constant temper- 
ature. The BD simulations were carried out for iVc=1113 
hard spheres, at fixed packing fraction c/)c=0.1, and for 
a number of droplets in the range iVd=10 to 667, i.e. at 
different droplets volume fractions (pd- To improve the 
statistical accuracy of the sampled quantities each state 
point was simulated four or eight times. 

The equilibrium phases of the model are investigated 
by carrying out simulations (for t = 1000 tb) at fixed 
(pc =0.1 and different values of (0d,e). The final phase 
for each state point is reported in Fig. [T] For both 
Cd/cc=0.75, and 1.5 we find a region where the fluid is 
stable (squares) and a region where gas-liquid phase sep- 
aration occurs (circles). We stress that this separation 
is not demixing of two species, but separation in a fluid 
that is rich in both droplets and colloids, and a dilute gas. 
Additionally, for crci/cc=1.5 we flnd a region where the 
fluid separates in a fluid and a binary crystal. In FiglTjD) 
a critical point and a triple point are recognizable. 

To study the cluster and gelation transitions wc choose 




Figure 2. (Color online) Simulation snapshots at i = 800 tb 
for (Td = 0.75. Droplets are shown as small red (dark grey) 
spheres while dilTerent clusters are colored (shaded) differ- 
ently, a) Fluid of clusters at $d=0.63% b) Fluid of clus- 
ters at $d = 1.05% c) Fluid at percolation point $d = 1.48% 
d)$d=1.9% Percolated state. 



state points within the gas-liquid envelope (see Fig. 8 of 
Ref. [i3)- The energy scale of the droplet-particle inter- 
action was fixed at e = IQQk'oT. Fig. [2] shows simulation 
snapshots at different droplets volume fractions (/)d, at 
t = 800 Tb . Different clusters are colored (shaded) differ- 
ently. A cluster is defined as a group of bonded particles. 
A bond between two particles is formed only when they 
share the surface of the same droplet. At low droplet 
volume fractions, as shown in Figs.[2^)-b), many clusters 
coexists with single particles and the system is therefore 
characterized by a high degree of dynamical heterogene- 
ity. At the higher droplet volume fractions (Figs. [2t)- 
d)) a large cluster spans the entire simulation box. A 
visual inspection of the snapshots indicates that a tran- 
sition from a fluid of clusters to a percolated network 
occurs, and that it is controlled by the amount of emul- 
sion droplets. 

The percolation transition was analyzed by carrying 
out a cluster analysis of eight independent runs at dif- 
ferent droplet volume fractions 0d. Figure [3^) shows the 
probability Pl that a particle belongs to the largest clus- 
ter as a function of the rescaled droplet volume fractions 
02 = </'d/</'d- The critical volume fraction c/i^ is defined 
as the volume fraction of the percolation transition, i.e. 
where Pl = 0.5. At small 0d ^^ observe a fluid of in- 
creasingly large clusters, until a sharp transition to a per- 
colated network is observed at large (j)^ . The percolation 
critical values are </)j=0. 01415, 0.039, for droplet sizes 
fd/o'c=0.75, and 1.5, respectively. The small amounts 
of secondary fluid necessary for the percolation transi- 
tion to occur are in good qualitative agreement with the 
results of Koos and Willenbacher [l| who found a transi- 
tion from a fluid to a gel behavior at a volume fraction of 
water (f> = 0.003 in a suspension of particles with volume 




Figure 3. (Color online) a) The probability Pl that a parti- 
cle belongs to the largest cluster as a function of the rescaled 
droplet volume fractions <^5 = 4'd/4'd- The critical values of 
the droplet volume fraction are 0^=0.01415, 0.039, for droplet 
sizes crd/crc=0.75, and 1.5, respectively, b) The radius of gy- 
ration of the clusters R^ versus the size n of the cluster, c) 
The time evolution of the number of clusters ncis(i)/ncia(0) 
for (j}*i — 0.75. The number of clusters at t = is ncis(0)=48. 



fraction (j) = 0.111. We find that the transition value 
depends on the relative size of droplets and particles. 
Therefore we expect an influence of the droplets poly- 
dispersity on the location of the percolation transition. 
More information about the percolation mechanism can 
be extracted from the analysis (Fig. IUd)) of the radius 

of gyration _Rg = ^\/X]j<7 ki ~ ^j\'^ o^ clusters of size 

n (number of constituents particles). For small sizes the 
clusters grow with a fractal dimension df ~ 1.80 ± 0.05; 
this value suggests three-dimensional diffusion limited 
cluster aggregation (DLCA) [18l - l20l |. At larger cluster 
sizes the aggregation is characterized by a fractal dimen- 
sion (if = 2.6±0.1,i.e. a random aggregation of different 
mesoscopic clusters, which form the different branches of 
the final percolating cluster. This value of d{ is also con- 
sistent with the fractal dimension di = 2.7 ±0.1 obtained 
from the slope of the cumulative integral of the radial 
distribution function of the percolated networks. Since, 
the crossover to random aggregation occurs at sizes that 
are about a quarter the lateral size of the simulation box, 
we do not expect finite size effects to dramatically change 
the location of the crossover. Lastly, in Fig. [3]:) the time 
evolution of the number of clusters ncis{t) is shown. Sur- 
prisingly, the number of clusters remains constant during 
the simulation time, indicating that the fluid of clusters 



is stable. 

We study dynamics via the colloid self-intermediate 
scattering function 



1 



Nr 



^^(*) = ^<E 



-'q-rj(t)„iqTj(0) 



(4) 



where Vj(t) is the position vector of particle j at time i, 
and q is a wave vector. The brackets indicate an aver- 
age over time origins and four independent simulations. 
The sampling was done for a time i/re = 6000. Fig. 0^) 
shows the behavior of Fq{t) as a function of t at fixed 
qac = 6.28 (yellow triangles) and qac = 0.6 (red circles) 
for different state points. For q<Tc = 6.28, (f)'^^0.75, and 
1.05, the function decays to zero. This behavior rep- 
resents the normal dynamics of a fluid, where particles 
quickly become uncorrelated (in time) at length scales of 
the order of a particle diameter. At 0^=1.35, i.e. well 
inside the percolation region, the decay is much slower. 
At qac — 0.6 (red circles) the transition from ergodic 
02=0.75 to non-ergodic behavior at 05=1.35 is a strong 
indication that a structural arrest characteristic of a gel 
transition [12| has occurred at percolation. Figure |3Jd) 
shows that all wave vectors remain ergodic at 0^=0.75. 
Therefore, the reason for the stability of the fluid of clus- 
ters (Fig. |3t)) is not slow dynamics, but the saturation 
of the attraction. The structural arrest at 0^=1.35 is 
clearly shown in Fig. Ht), with a plateau (non-ergodic 
behavior) that develops fully only for small qa^- 

The mechanical response of the system is characterized 
via the shear stress autocorrelation function Cs(t^ 

cs{t) = -^(5„^(0)5„^(i)) (5) 

where V is the system volume and the brackets indicate 
an average over time origins and eight independent sim- 
ulations. The sampling was done for I/tb = 400 and 
t/TB = 1200. The function Saj3{t) represents the off- 
diagonal elements of the stress tensor 
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where r^ is the a-component of the distance vector r^j 
between particles i, and j, and fi^ is the /3-component of 
the pairwise force. Fig. [5^) shows the stress autocorrela- 
tion function as a function of time. In a fluid, (0^ = 0.03) 
Cs (t) quickly decays to zero with a stretched exponential 
form. The transition to a power law decay (linear region 
in the log-log plot for 0^ = 0.45, 0.75) is a sign of an 
approaching gel transition lll|. When a percolated net- 
work is formed (0^ = 1.05, 1.35, 1.5) a plateau in the cor- 
relation function develops. The plateau is an indication 
that the system is able to sustain stress. At equilibrium, 
the time integral of c^(t) gives the fluid viscosity via the 
Green-Kubo relation uM 
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Figure 4. (Color online) Self-intermediate scattering function 
F[q, t) as a function of time t. a) For guc = 6.28 (yellow 
curves) and qoc = 0.6 (red curves) for scaled droplet volume 
fractions ^J^— 0.75, 1.05, 1.35. b) For (/)5=0.75, and go-c=0.6, 
1.2, 1.8, 3.8, 4.2, and 8.0, from top to bottom, c) Same as b), 
but for 05 = 1.35. 
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In Fig. [5]d) the values of r}{ty^ax) are plotted against the 
scaled droplet volume fraction (p*^ for two different droplet 
sizes crd/crc=0.75, and 1.5. The transition from liquid- 
like to solid-like behavior is not sharp because of the 
finite time used in the simulation. We stress that the 
values plotted in Fig. [Sb) represent the viscosity only if 
the stress correlation function decays to zero. Nonethe- 
less, the value of the integral ([7]) is an indicator of the 
changing mechanical response of the system and it shows 
that small droplets influence the flow behavior more than 
large droplets. 

In this article a model was introduced that describes 
the formation of a fluid of stable clusters and gels when 
emulsion droplets are added to a hard spheres suspension. 
The particle aggregation into clusters is due to the Pick- 
ering effect, while the mechanism for cluster-cluster ag- 
gregation is due to colloids bonding to multiple droplets. 
The cluster-cluster aggregation leads to the formation of 
a spanning network of colloidal particles via a percolation 
transition. A crossover from DLCA to random aggrega- 
tion of large clusters is observed, and the critical percola- 
tion density of emulsion droplets depends on the droplet 
to particle diameter ratio. The analysis of the dynamics 
via the colloid self-intermediate scattering function shows 
a strong structural arrest at small wavelengths, indicat- 
ing the occurrence of a gel transition in correspondence 
of the percolation transition. The analysis of the stress 
correlation function also suggests the ability of the per- 
colated networks to sustain stress. Therefore, both the 
dynamical and mechanical properties of these networks 
arc those of gels, in agreement with the experiments [l[. 
Due to bond saturation there is a lack of available at- 
traction in the final configurations that, combined with 
a strong particles-droplet bonding, suggests that aggre- 
gation and collapse of clusters and gels is hindered even 
at very long timescales. 

Furthermore, we traced the stability phase diagram for 
the model and found that both gas-liquid and fluid-solid 
separations are possible. In particular, we found that 
binary crystals can self-assemble at high concentration of 
emulsion droplets. These phases were not yet observed 
experimentally in Pickering emulsions. Although their 
realization represent a big technical challenge, emulsions 
could prove to be a novel way to control the formation of 
complex colloidal crystal structures. 



Figure 5. (Color online) a) Stress autocorrelation func- 
tion Cs(i)/cs(0) at (Td/o"c=0.75 and different values of 
scaled droplet volume fraction (j)*^ b) Scaled viscosity fj = 
ri{tniax)/'ri{'tinax,4>d = 0.03) for two different droplet sizes 
o"d/o"c=0.75, and 1.5. 
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